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The effect of dephasing on electron transport through a benzene molecule is carefully examined 
using a phenomenological model introduced by Biittiker. Within a tight-binding framework all the 
calculations are performed based on the Green's function formalism. We investigate the influence of 
dephasing on transmission probability and current- voltage characteristics for three different configu- 
rations (ortho, meta and para) of the molecular system depending on the locations of two contacting 
leads. The presence of dephasing provides a significant change in the spectral properties of the 
molecule and exhibits several interesting patterns that have so far remain unexplored. 
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I. INTRODUCTION 

Ongoing trend of miniaturizing electronic devices even- 
tually approaches the ultimate limit where even a single 
molecule can be used as an electrical circuit element. Idea 
of devicing a single molecule as the building block of fu- 
ture generation electronics seems fascinating because of 
the possibility to assemble a large number of molecules 
onto a chip i.e., remarkable enhancement in integration 
density can take place Discovery of sophisticated 
molecular scale measurement methodologies such as 
scanning tunneling microscopy (STM), atomic force mi- 
croscopy (AFM), scanning electro-chemical microscopy 
(SECM), etc. have made it possible to study electron 
transport phenomena in molecular bridge systems 0. 
Unlike conventional semiconductor electronic transport 
properties can not be investigated with Boltzmann trans- 
port equation as quantum coherence has a crucial signifi- 
cance on electron transport at this length scale and a full 
quantum mechanical treatment is needed [J 3 ■ At low 
temperature and small bias electron transport is quite 
successfully explained within the framework of Landaucr 
formalism [||. But, validity of Landauer formulation is 
limited within coherent transport regime where it is as- 
sumed that electrons suffer only elastic scattering within 
the conductor. This assumption allows us to neglect any 
kind of inelastic processes, which is rather unrealistic un- 
less the time scale related to transport phenomena is 
much faster than the nuclear motion . 

A main source of dephasing is electron-phonon (e-ph) 
interaction in molecular transport junctions. Experi- 
mentally the strength of e-ph interaction in molecular 
systems can be quantified through the measurement of 
molecular vibrational spectrum using inelastic tunnel- 
ing spectroscopy @-H[- Roughly, it corresponds to find- 
ing locations of the peaks in the second order derivative 



of a current-voltage characteristic, where the voltages 
at the peaks match with eigenenergies of the phonons. 
Rigorous analysis of e-ph interaction can be done self- 
consistcntly with the help of density functional theory 
(DFT) and non-equilibrium Green's function formalism 
(NEGF) [TO, [ll|. But this method is notably time con- 
suming and becomes very difficult to do even for a molec- 
ular conductor comprising a moderate number of atoms. 
Another way out was proposed by Biittiker over two 
decades ago [HI, • He came up with an elegant idea of 
incorporating the effect of phase breaking due to inelas- 
tic scattering phcnomenologically by introducing some 
fictitious voltage probes into the coherent system. While 
moving from source to drain electrons get scattered into 
the fictitious floating voltage probes and then re-emitted 
into the device by loosing the phase memory. To im- 
plement this idea each inelastic scatterer is modeled in 
terms of electron reservoir (or infinite impedance volt- 
age probe) and coupled to the device. Since the net 
current flowing through the fictitious probes has to be 
zero to satisfy the conservation of total number of par- 
ticles, electro-chemical potentials (fi) of each lead has to 
be adjusted accordingly. Due to its appealing simplic- 
ity Biittiker probe model has been extensively used in 
studying quantum transport in low-dimensional systems. 

The idea of using molecules as active components of a 
device was s ugg ested by Aviram and Ratner over three 
decades ago [14|. Since then several ab-initio and model 
calculations have been performed to investigate molecu- 
lar transport theoretically floT - 1231 ]. But experimental re- 
alizations took a little longer time to get feasible. In 1997, 
Reed and co-workers studied the current-voltage (I-V) 
characteristics of a single benzene molecule attached to 
electrodes via thiol groups [24| . Various other experi- 
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ments have been reported in literature exploring many 
interesting features e.g., ballistic transport, quantized 
conductance, negative differential resistance (NDR), gate 
controlled transistor operation to name a few. Effect of 
dephasing in molecular systems has been studied in few 
literatures both theoretically [Ifl [2(| as well as cxperi- 
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mentally. In 1989, Amato and Pastawski investigated the 
effect of inelastic scattering processes on a tight-binding 
linear disordered chain using Buttiker probe model . In 
2007, Datta et al. proposed a different phenomenological 
model which provides the flexibility of adjusting degree 
of phase and momentum relaxation independently [271 ]. 
In 2008, Nozaki et al. obtained transport properties 
through various molecular junctions by Buttiker model 
using Extended Hiickel theory (28|. In 2009, Guo and 
co-workers have reported an ab-initio calculation (29j in- 
cluding the effect of phase breaking in 1,4-benzenedithiol 
molecule attached to two Al electrodes, gold QPC and 
a very thin Al(001) nanowire. They have also done a 
comparative analysis between Buttiker probe model and 
Datta dephasing model [13] ■ Though various efforts have 
already been made to explain the basic features of elec- 
tron dephasing on molecular transport, but to the best 
of our knowledge, no rigorous effort has been made so far 
to unravel the combined effect of quantum interference 
and electron dephasing on molecular transport. This is 
the main motivation behind this work. 

In our present article we study two terminal elec- 
tron transport through a single benzene molecule includ- 
ing the effect of dephasing by a discrete tight-binding 
model. We obtain the transmission probability us- 
ing Green's function technique within the framework of 
Landaucr-Buttikcr formulation. Transmittance-energy 
and current-voltage characteristics are obtained for three 
different configurations (para, ortho and metd) of ben- 
zene molecule depending on three different positions of 
the drain with respect to the source. The presence of 
dephasing provides a significant change in the spectral 
properties of the molecule and here we essentially focus 
our results in this aspect. Our model calculation can be 
extended further for any complicated molecular struc- 
ture. 

In what follows, we present the results. In section II, 
the model and the theoretical formulation are presented. 
Section III describes the results, and finally, we make our 
conclusions in section IV. 



an infinite impedance voltage probe (shown by the pink 
shaded area in Fig. [I} via perfect leads. The key idea is 
that the fictitious floating probes extract electrons from 
the device to the reservoir, but as the net current flowing 
through these probes is zero, so they re-inject the elec- 
trons into the conductor after phase randomization and 
thus effectively play the role of phase breaking scatterers. 

The Hamiltonian for the entire system is given by the 
sum of four terms 

H = H mo i + Hi ea d s + if tun + Hdephase- (1) 

The first term represents the Hamiltonian for the single 
benzene molecule consisting of six sites (N = 6) which 
is coupled to two electron reservoirs through conducting 




Benzene Molecule 

FIG. 1: (Color online). Schematic diagram of a single benzene 
molecule attached to two leads, namely, source and drain. 
Each molecular site is assumed to be coupled to a fictitious 
electron reservoir, shown by the pink shaded area, which takes 
into account the effect of phase breaking processes. 

leads, namely, source and drain. The molecule is modeled 
by the nearest-neighbor tight-binding Hamiltonian which 
in Wannicr basis can be written as, 

Hrnol = 6d l d i + V i4+A + h.C.], (2) 

i i 



II. MOLECULAR MODEL AND THEORETICAL 
FORMULATION 

We start by describing our model, illustrated in Fig. [TJ 
where a single benzene molecule is connected symmet- 
rically to two one-dimensional, semi-infinite leads, com- 
monly known as source and drain. The leads are char- 
acterized by the electrochemical potentials us and \xb, 
respectively, under the non-equilibrium condition when 
an external bias voltage is applied. Both the two leads 
and the molecule are simulated by a simple tight-binding 
Hamiltonian within nearest-neighbor hopping approxi- 
mation. Following Biittiker's idea, the effect of dephasing 
due to inelastic scattering is incorporated in the model 
phenomcnologically by coupling each spatial sites to a 
different fictitious electron reservoir or in other words 



where e refers to the site energy of an electron at the 
i-th site of the molecular system and v represents the 
isotropic nearest-neighbor coupling strength between the 
molecular sites. d\ and di correspond to the creation and 
annihilation operators, respectively, of an electron at the 
i-th site of the molecule. 

Similarly the second and third terms of Eq. ((TJ) denote 
the Hamiltonians for the one-dimensional semi-infinite 
leads (source and drain) and molecule-to-lead coupling. 
In Wannier basis representation they can be written as 
follows. 

Hieads = Hs + Hp 

= \ J2 6lC n Cn + t^li+l " + h - C '} \ ' 

a=S,D In n ) 

(3) 
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and, 



H t un — Hs,mol + He, 

.mot 

= T S [d\c Q + h.c] + r D [dlc N+1 + h.c.]. (4) 

Here, ei and v stand for the site energy and nearest- 
neighbor hopping between the sites of the two leads. c? n 
and c ra are the creation and annihilation operators, re- 
spectively, of an electron at the site n of the leads. The 
hopping integral between the source and the molecule is 
rs, while it is td between the molecule and the drain. 
Usually benzene molecules are attached to gold leads via 
thiol groups (-SH bond) in the chemisorption technique. 
Here, the effect of such substituent (e.g., thiol group) is 
included in the parameters ts and td- The source is as- 
sumed to be connected at the site number 1, while the 
drain is coupled to the p-th site, which is variable, of the 
molecular system. The integer p becomes 2, 3 and 4 for 
the ortho, meta and para configurations, respectively 

In order to incorporate the effect of dephasing, exter- 
nal electron reservoirs arc phcnomenologically introduced 
and they arc connected to all the molecular sites (six) 
through fictitious probes. The phase breaking probes 
are identically represented within tight-binding frame- 
work like source and drain as given below. 



H, 



dephase 



^ ^ probes jjd,tun 



i = l 



(5) 



where, 



Tjd,probes \ "> Dii Di , \ * , r Dii Di , > l t a \ 

m m 

and, 



TT d,tun Dii » . j 

H i ' =r] c ] d l + h.c. . (7) 



Here, td and td being the site energy and the hopping 
strength between the nearest-neighbor sites of the de- 
phasing probes, respectively, c^(c^) is the creation 
(annihilation) operator of an electron at the m-th site of 
the i-th lead. The dephasing strength is characterized by 
the coupling parameter 77. 

To obtain the transmission probability of an electron 
through such a molecular bridge system, we use Green's 
function formalism. Within the regime of coherent trans- 
port and in the absence of Coulomb interaction this tech- 
nique is well applied. 

The single particle Green's function operator repre- 
senting the entire system for an electron with energy E 
is defined as, 



G={E-H + iri) 1 



(8) 



where, rj — > + . 

Following the matrix form of H and G the problem 
of finding G in the full Hilbert space H can be mapped 



exactly to a Green's function G^f/j corresponding to an 
effective Hamiltonian in the reduced Hilbert space of the 
molecule itself and we have, 



Q=G e JL = 



N 

E — H mo i — Eg — Tid — Sfc 

fc=i 



where, 
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S,mol(mol,R)G S(D)H S,mol{mol,R) 
T_rd,tun1[ Trd,tun 



(9) 



(10) 



These Ss(r>) an d E& are the contact self-energies intro- 
duced to incorporate the effect of coupling of the molecule 
to the source (drain) and the fictitious electron reservoirs. 
It is evident from Eq. (jTUJ) that the form of the self- 
energies are independent of the molecule itself through 
which transmission is studied. Using Dyson equation the 
analytic form of the self energies can be evaluated as fol- 
lows, 



(11) 



E — e v - ov 



Here, a v = (E - e„)/2 - i^/t 2 v - (E-e„) 2 /4; and v = 
5, D, 1, 2, . . . N. For source and drain, v = S{D), r„ = 
T S(D), e v = e, and i„ — t; while considering the dephasing 
probes, v = k (k = 1,2,... N), t v = rj, t v = e;, and 
t v = td- 

Following D'Amato-Pastawski model [a] the effective 
transmission probability of an electron from source to 
drain including the partially phase breaking process is 
given by the following expression. 



T, 



eff 



T S .D + 



N 

E' 



w 



• ■J 



■3,S- 



(12) 



The above expression has a significant physical impli- 
cation. The first term is the coherent contribution to 
electron transmission, whereas the second term refers to 
the incoherent component of tunneling due to electrons 
suffering phase randomizing scattering at the fictitious 
dephasing reservoirs. Ti j is the transmission probability 
between any pair of reservoirs i, j = S, D, 1, 2 . . . N) 
and it is expressed as follows, 



T iij =Tr[T i g r T j g a ]. 



(13) 



Ti's are the coupling matrices representing the coupling 
between the molecule and the leads and they are mathe- 
matically defined by the relation, 



r, = i [sj - ej] 



(14) 



Here, E£ and E£ are the retarded and advanced self- 
energies associated with the fc-th lead, respectively. 

In the second term of Eq. (O, W matrix is defined 
by the following relation, 



^• = [(1-^,0^-^(1-^)], 



(15) 
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where, 



(16) 



Ri i is the reflection probability of an electron from i-th 
lead. 

It is shown in literature by Datta et al. 0, Q that the 
self-energy can be expressed as a linear combination of a 
real and an imaginary part in the form, 



(17) 



The real part of self-energy describes the shift of the en- 
ergy levels and the imaginary part corresponds to the 
broadening of the levels. The finite imaginary part ap- 
pears due to incorporation of the semi-infinite leads hav- 
ing continuous energy spectrum. Therefore, the coupling 
matrices can easily be obtained from the self-energy ex- 
pression and is expressed as, 



-2Im(Ei). 



(18) 



Considering linear transport regime, at absolute zero 
temperature the linear conductance (<? e //) is obtained 
using two-terminal Landaucr conductance formula, 



9eff 



2e 2 

— T eff{ E F)- 



(19) 



With the knowledge of the effective transmission proba- 
bility we compute the current-voltage (I-V) characteris- 
tics by the standard formalism based on quantum scat- 
tering theory. 



T eff [f s (E)-f D (E)]dE. (20) 



is the Fermi func- 



Here, f s{D )(E) = 14 

tion corresponding to the source and drain. At absolute 
zero temperature the above equation boils down to the 
following expression. 



TeffdE. 



(21) 



sandwiched between source and drain. The orbital ener- 
gies (e) of the molecule are set to zero and the nearest- 
neighbor hopping matrix element (t) is fixed at 2.4 eV [TtJ 
for an aromatic ring structure of the benzene molecule. 
The other hopping parameters in the leads (t) and the 
dephasing probes (td) are fixed at 4eV, while the site en- 
ergies arc chosen to be zero i.e., e/ = = 0. The dephas- 
ing strength, characterized by the parameter 77, is fixed 
at 1 eV. Throughout our analysis the coupling strengths 
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FIG. 2: Variation of transmission probability as a function of 
energy for a benzene molecule where the drain is connected 
at the para position with respect to the source. The solid and 
dashed curves correspond to the results in the absence and 
presence of dephasing, respectively. 

of the molecule to the source and drain, characterized by 
the parameters t$ and 775 , are set at 1 eV. The molecular 
coupling to the contacting leads is one of the most im- 
portant factors which regulates the electronic transmis- 
sion through a molecular wire and this phenomenon has 
been extensively studied in our previous works [22|, [23| . 
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FIG. 3: The results of Fig. [2] are re-plotted in a logarithmic 
scale to clarify the effect of dephasing more clearly. 



In our present work we use the above expression assuming 
that the potential drop takes place only at the boundary 
of the molecular system. This assumption is good enough 
for molecules of smaller size. Throughout our study we 
choose c = e = h = 1 for the sake of simplicity. 



III. NUMERICAL RESULTS AND DISCUSSION 

In this section we present the results obtained by nu- 
merical simulation considering a single benzene molecule 



Accordingly, here we describe all the essential features 
of electron transport for a particular molecular coupling 
strength. The molecular transport properties are also 
highly sensitive on the molecule-to-lead interface geome- 
try, and, in the present work we illustrate how they are in- 
fluenced in the presence of electron dephasing. To reveal 
these facts first we discuss the behavior of transmittance- 
energy characteristics and then we describe the nature of 
current- voltage spectra for three different configurations 
(ortho, meta and para) of the benzene molecule. 

Here, we establish all the results at absolute zero tern- 
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perature and restrict ourselves within elastic dephasing 
so that no energy is exchanged between the transported 
electrons and the external fictitious reservoirs, only the 
phase informations get lost. This condition allows the 
incoherent electrons to generate a steady-state current 
through the sample. 

A. Transmittance-energy characteristics 

1. Molecule coupled symmetrically 

In Fig. [5] we show the variation of transmission prob- 
ability as a function of injecting electron energy for a 
benzene molecule, where the molecule is coupled sym- 
metrically to the source and drain i.e., the upper and 



cific energies, whereas it almost vanishes for all other en- 
ergy values. At these resonances, the transmission prob- 
ability T e ff goes to unity. All the resonant peaks are 
associated with the energy eigenvalues of the benzene 
molecule, and therefore, we can say the transmittancc 
spectrum is a fingerprint of the electronic structure of the 
molecule. The situation becomes much interesting when 
the effect of dephasing is introduced. It shows that the 
magnitude of the resonant peaks gets suppressed enor- 
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FIG. 4: Transmission probability versus energy characteris- 
tics of a benzene molecule where the drain is connected at 
the ortho position with respect to the source. The solid and 
dashed lines represent the results in the absence and presence 
of dephasing, respectively. 

lower arms of the molecular ring have identical length. 
This is the so-called para configuration. The solid curve 
represents the variation of coherent transmission proba- 
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FIG. 5: Transmission probability versus energy character- 
istics of a benzene molecule where the drain is coupled to 
the meta position with respect to the source. The solid and 
dashed curves correspond to the results in the absence and 
presence of dephasing, respectively. 

bility, while the dashed curve depicts the result in the 
presence of electron dephasing. It is observed that, in 
the absence of dephasing the transmission probability ex- 
hibits sharp resonant peaks (see solid line) for some spe- 
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FIG. 6: The results of Fig. [4] are re-plotted in a 'log' scale. 

mously (dashed line) compared to the coherent trans- 
mission (solid line). This is due to the increased rate 
of scattering in the additional fictitious electron reser- 
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FIG. 7: The results of Fig. [5] are re-plotted in a 'log' scale. 

voirs. Additionally, we also notice that the widths of the 
resonant peaks become broadened. This broadening is 
caused by the dominance of phase randomizing effect over 
backscattering due to loss of phase coherence of interfer- 
ing electrons in the presence of dephasing reservoirs. To 
judge the effect of dephasing more transparently, in Fig. [3] 
we re-plot the results of Fig. |5| in a logarithmic scale, 
where the solid and dashed lines represent the identical 
meaning as above. Since the transmission probability lies 
within the range to 1, it becomes zero or negative in the 
logarithmic scale as shown in Fig. [3] It is observed that 
the variations in the transmission peaks between the two 
curves (solid and dashed) get suppressed, but the broad- 
ening effect in the presence of dephasing can be much 
clearly observed from the 'log' scale. 
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2. Molecule coupled asymmetrically 

A significant change in the transmission spectrum is 
observed when the benzene molecule is coupled asym- 
metrically to the source and drain i.e., the upper and 
lower arms of the molecular ring have unequal lengths. 
For the ortho configuration the results are given in Fig. 0] 
while in Fig. [5] the results are shown for the meta con- 
figuration. From the spectra (Figs. 0] and [S]) we notice 
that in the absence of dcphasing some resonant peaks 
do not reach to unity and their amplitudes are signif- 
icantly reduced compared to the other resonant peaks. 
This is solely due to the effect of quantum interference 
between the electronic waves passing through the upper 
and lower arms of the molecular ring. Quite interestingly 
we observe that in these asymmetric molecular wires {or- 
tho and meta) two anti-resonant states appear in the 
transmittance spectrum where transmission probability 
drops exactly to zero. These anti-resonant states are spe- 
cific to the interferometric nature of the molecular sys- 



transmission peaks in presence of dcphasing is much pro- 
nounced in the case where the drain is connected at the 
meta position of the benzene molecule rather than the 
ortho position, and, it can be clearly observed from the 
Log e T e ff vs. E plot instead of the transmittance-energy 
spectra. 



B. Current-voltage characteristics 

All the above features of electron transmission become 
much clear from our study of current-voltage (I-V) char- 
acteristics. The current through the molecular system 
is obtained by integrating the transmission function fol- 
lowing Eq. (f2Tj) . Here we fix the Fermi energy Ep at 
zero. 



1. Molecule coupled symmetrically 
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FIG. 8: (Color online). I-V characteristics of a benzene 
molecule where the drain is connected at the para position 
with respect to the source. The orange and green curves 
correspond to the currents in the absence and presence of 
dephasing, respectively. 

tern and their positions in the energy scale are indepen- 
dent of the molecule-to-lead coupling strength. A careful 
observation reveals that a sharp peak is followed by a 
sharp dip and vice versa across the anti-resonant energy 
which is quite analogous to a Fano-like line shape. All 
these anti-resonant states disappear as long as dephas- 
ing of electrons is taken into account those are clearly 
observed from the dashed lines of Figs. 0] and [SJ Here 
also the magnitudes of the resonant peaks get decreased 
and they become broadened as we describe earlier in the 
para configuration. The results of Figs. |4] and [5] are re- 
drawn in Figs. [5] and respectively, where the transmis- 
sion probabilities are determined in the 'log' scale. From 
these spectra (Figs. [6] and [7|) the characteristic features of 
anti-resonant states in the case of coherent transmission 
and the broadening of resonant peaks in the presence of 
electron dephasing are clearly noticed. Our results also 
predict that in the asymmetric molecular ring the po- 
sitions of the anti-resonant states are strongly sensitive 
on the location of measuring electrodes. Broadening of 



In Fig. |8] we plot the variation of current as a func- 
tion of applied bias voltage for a benzene molecule where 
the leads are connected symmetrically. The orange and 
green lines represent the currents in the absence and 
presence of electron dephasing. The current varies quite 
continuously with the applied bias voltage V . Depending 
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FIG. 9: (Color online). I-V characteristics of a benzene 
molecule where the drain is connected at the ortho position 
with respect to the source. The pink and deep blue lines rep- 
resent the currents in the absence and presence of dephasing, 
respectively. 

on the molecular coupling to the side attached leads, the 
current shows continuous- like or step-like behavior [22| . 
Tuning the molecular coupling strength the current am- 
plitude through the molecular wire can be regulated 
nicely for a fixed bias voltage. This provides an inter- 
esting behavior in designing molecular electronic devices. 
Both in the absence and presence of electron dephasing 
current shows almost identical variation though a signif- 
icant change is observed in their conductance spectra. 
This is due to the fact that the dephasing broadens the 
transmission peaks, while it also suppresses the magni- 
tude. These two effects nullify each other showing almost 
identical spectrum in the I-V characteristics. 
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2. Molecule coupled asymmetrically 

The effect of dephasing becomes much prominent in 
the current-voltage characteristics when the molecule is 
coupled asymmetrically to the source and drain. For the 
ortho configuration, the currents are shown in Fig. |H1 
while in Fig. [10] the results are shown when the molecule 
is coupled to the leads in the meta configuration. Quite 
similar to the para configuration (Fig. [8]), in these two 
cases {ortho and meta configurations) also the current 
varies almost continuously with the applied bias voltage 
V. In these asymmetric configurations, the current am- 
plitude gets enhanced significantly when the phase break- 
ing effect is considered (deep blue curve in Fig. [9] and 
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FIG. 10: (Color online). I-V characteristics of a benzene 
molecule where the drain is connected at the meta position 
with respect to the source. The light blue and red lines corre- 
spond to the currents in the absence and presence of dephas- 
ing, respectively. 

light blue curve in Fig. ITU)) compared to the case where 
the molecule is free from any phase randomizing process 
(pink line in Fig. |9] and red line in Fig. [TO]) . Such an 
enhancement in current amplitude takes place because 
of notable broadening of the transmission peaks due to 
reduced destructive interference for dominance of phase 
decoherence over the backscattering effect. 



IV. CLOSING REMARKS 

To summarize, in the present work we have examined 
electron transport phenomena through a single benzene 
molecule including the effect of phase breaking using a 
phenomenological model introduced by Buttiker. Exper- 
imentally the benzene molecule is usually coupled to no- 
ble metal electrodes (e.g., gold leads) through thiol link- 
ing groups. In our theoretical work, we present the en- 
tire system (source-molecule-drain) by a discrete lattice 
model within a simple one-electron, tight-binding frame- 
work. We have studied transmission-energy spectrum 
and current- voltage characteristics considering three dif- 
ferent configurations of the system depending on the po- 
sition of the drain contact (ortho, meta and para posi- 
tions). We have observed that the sharp transmission 
peaks are broadened and their magnitude is also sup- 



pressed in the presence of dephasing due to loss of phase 
coherence and thus it provides an enhancement in the 
current amplitude. But the broadening effect and so as 
the increase in current amplitude is much pronounced in 
the case where the drain is connected asymmetrically to 
the molecule with respect to the source. Another sig- 
nificant feature is that the two anti-resonant states, ob- 
served in the transmittance spectrum of an asymmetri- 
cally connected benzene molecule, completely disappear 
in the presence of strong dephasing due to reduced de- 
structive interference. Thus the effect of dephasing can 
provide a significant change in the current-voltage rela- 
tion. 

In the Buttiker probe model fictitious voltage probes 
are introduced into the coherent system. These virtual 
probes introduce additional resistance in the system by 
relaxing the momentum of the electrons by partially re- 
flecting them i.e., back-scattering. This model is appro- 
priate for those molecular systems where the dominant 
scattering source does not conserve electron momentum, 
for example, electron-phonon scattering which is very 
important in molecular conductor in finite temperature. 
However, Buttiker model does not consider the effect of 
momentum relaxation. In experimental situations de- 
phasing may also arise from electron-electron scattering 
which destroy phase but not momentum [27|, H^l ■ Datta 
model presents the flexibility of tuning phase and mo- 
mentum independently of each other and we intend to 
study this phenomenological model in near future. 

Recent experimental progresses have significantly in- 
spired in model calculations and numerical simulations 
to study vibronic effects on electron transport in molec- 
ular junctions [30Tj33j . Most of the studies of vibronic 
effects in the resonant tunneling regime (for higher volt- 
ages) have illustrated that vibrational motion may affect 
the current-voltage spectra quite significantly. But, the 
point is that in the resonant tunneling region vibronic 
effect becomes significant only when the molecule is cou- 
pled to the measuring leads weakly. In a very recent 
work [34[ Benesch et al. have studied vibronic effects on 
conductance in molecular junctions where a benzene ring 
is coupled weakly to two gold electrodes and shown that 
the vibronic motion affects current- voltage spectra. On 
the other hand, when the molecule is coupled strongly to 
the electrodes, the vibrational effects in resonant trans- 
port are almost negligible. Thisphenomenon has been 
clearly justified in the reference [351 ] . In the present work 
we have investigated all the essential features of electron 
transport in the limit of strong molecule-to-electrode cou- 
pling, and accordingly, we neglect the vibronic effect. 

Throughout our work we have neglected the intra- and 
inter-site Coulomb interactions as well as the effect of 
the leads, which we wish to consider in our future stud- 
ies. Another important assumption is the zero temper- 
ature approximation. Although the electron-phonon in- 
teraction becomes significant at non-zero temperature, 
but in our case, we have already considered the dephas- 
ing effect phcnomcnologically at every molecular sites, 
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which takes into account the effect of dcphasing of elec- 
trons without any exchange of energy between the elec- 
tron and the external leads. Accordingly, we compute 
the current-voltage relation at absolute zero tempera- 
ture. At finite non-zero temperature, the transmission 
peaks will broaden due to thermal broadening effect, 
but the basic features will not change significantly as 
long as the thermal energy (ksT) is less than the av- 
erage level spacing of the benzene molecule. Now-a-days 
various organic compounds and polymers (for example, 



Tris(8-hydroxyquinolino) aluminum (Ak}3) and 1,2,3,4,5- 
pentaphenylcyclopentadiene (PPCP)) are used as elec- 
troluminescent devices. Knowledge of electronic struc- 
ture and current-voltage relation is imp ortant to eluci- 
date the light emitting mechanism [361 ] . It may be in- 
teresting to see the effect of electron dephasing on this 
mechanism. Our presented results may be useful in un- 
derstanding the effect of phase breaking processes on the 
basic two-terminal molecular transport phenomena. 
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